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Abstract 


Despite  very  similar  protein  structures,  ascorbate  peroxidase  (APX)  and  yeast  c)rtochrome-c- 
peroxidase  (CCP)  stabilize  different  radical  species  during  enzyme  turnover.  Both  enzymes 
contain  similar  active  site  residues,  including  the  tryptophan  that  is  oxidized  to  a  stable  cation 
radical  in  CCP.  However,  the  analogous  trytophan  is  not  oxidized  in  APX,  and  the  second 
oxidizing  equivalent  is  retained  as  a  porphyrin  tt-cation  radical.  In  this  study,  we  provide  an 
improved  computational  approach  to  estimate  the  contribution  of  solvent  and  protein 
electrostatics  to  the  energetics  of  tryptophan  cation  radical  formation  in  the  two  enzyme 
environments.  The  Protein  Dipoles  Langevin  Dipoles  (PDLD)  model  is  combined  with 
molecular  dynamics  to  estimate  the  role  of  discrete  solvation,  atomic  polarizabilities,  and 
dynamic  motional  averaging  on  the  electrostatic  potentials.  The  PDLD  model  shows  that  the 
protein  environment  of  CCP  stabilizes  the  tryptophan  cation  radical  by  330  mV  relative  to  that 
in  APX.  Analysis  of  the  components  contributing  to  this  difference  supports  proposals  that  the 
cation  binding  site  contributes  to,  but  is  not  the  sole  cause  of,  the  different  sites  of  radical 
stabilization.  The  enzymes  have  thus  evolved  this  distinction  using  several  contributing 
interactions  including  the  cation  binding  site,  solvent  access,  and  subtle  differences  in  protein 
structure  and  dynamics. 
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1.  Introduction 


Electron  transfer  occupies  a  position  of  central  importance  in  biological  systems.  Understanding 
biological  electron  transfer  requires  an  understanding  of  the  proteins  and  enzymes  with  redox  active 
cofactors  precisely  tuned  to  optimize  (both  thermodynamically  and  kinetically)  their  function  [1]. 
The  electrostatic  environment  of  such  a  cofactor  can  extend  significant  control  over  its  properties. 
If  the  cofactor  incorporation  into  two  proteins  is  similar  (e.g.,  the  same  ligand  set  for  a  metal  center), 
this  control  of  redox  functions  by  the  protein  should  be  amenable  to  computational  modeling  that 
includes  cofactor,  protein,  and  solvent  elements. 

The  peroxidases  comprise  a  class  of  enzymes  that  catalyze  the  H2O2  mediated  oxidation  of  an 
enormous  range  of  biological  substrates  [2,  3].  The  heme  containing  enzyme 
cytochrome-c-peroxidase  (CCP)  from  yeast  reacts  with  H2O2  to  form  a  two-electron  oxidized 
intermediate:  compound  ES.  This  intermediate,  in  turn,  accepts  electrons  sequentially  from  two 
ferrous  cytochromes  c.  One  oxidizing  equivalent  in  compound  ES  is  stored  as  an  oxyferryl 
(Fe(IV)=0)  heme,  while  the  other  is  stored  as  a  radical  species  on  an  amino  acid  side  chain.  In  CCP, 
tryptophan-191  is  preferentially  oxidized  to  a  Tc-cation  radical  to  house  the  second  oxidizing 
equivalent  [4-9].  Tryptophan- 191  in  CCP  is  positioned  as  part  of  a  so-called 
tryptophan-aspartate-histidine  triad  (the  histidine  being  the  heme  ligand)  thought  to  contribute  to  the 
stability  of  the  indole  cation  radical  on  tryptophan-191  [10].  Ascorbate  peroxidase  (APX),  a  plant 
peroxidase  that  has  a  very  similar  molecular  structure  to  CCP,  including  the 
tryptophan-aspartate-histidine  triad,  nevertheless  does  not  exhibit  an  indole  radical  upon  reaction 
with  hydrogen  peroxide  [11-13].  Listead,  a  porphyrin  ir-cation  radical  is  formed  [14].  This 
difference  in  location  of  the  stable  second  oxidation  equivalent  may  be  the  cause  of  the  different 
substrate  preferences  of  the  two  peroxidases:  large  proteins  for  CCP  (e.g.,  cytochrome  c)  and  small 
molecules  for  APX  (e.g.,  ascorbate).  It  is  believed  that  in  many  other  plant  peroxidases,  this  is 
accomplished  by  replacing  the  triad  tryptophan  with  another  residue,  such  as  phenylalanine,  which 
is  not  so  easily  oxidized.  In  this  report,  we  model  the  oxidation  of  tryptophan  using  the  x-ray  crystal 
structures  of  CCP  and  APX,  atomic  partial  charges  determined  ab  irutio,  and  an 
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electrostatic/molecular  dynamics  program  to  model  the  protein  and  solvent.  Calculations  of  the 
redox  potential  of  the  indole-to-indole  cation  radical  oxidation  in  tryptophan- 191  of  CCP  and  in  the 
analogous  tryptohan- 1 79  of  APX  reveal  that  the  redox  potential  of  the  APX  indole  oxidation  is  ~  330 
mV  higher  than  the  equivalent  reaction  in  CCP.  This  is  sufficiently  higher  in  APX  to  effectively 
account  for  the  absence  of  an  observed  indole  cation  radical  in  APX.  A  bound  monovalent  cation 
observed  in  the  APX  structure  ~9A  from  the  triad  indole  and  absent  in  the  CCP  structure  is  partly, 
though  not  completely,  responsible  for  the  calculated  increase  in  indole  redox  potential  in  APX 
relative  to  CCP,  with  other  more  stmcturally  subtle  factors  comprising  the  remainder  of  the 
difference.  These  calculations  build  on  earlier  modeling  studies  [11-13],  which  used  continuum 
models,  which  came  to  the  same  basic  conclusions:  There  is  an  electrostatic  origin  of  the 
destabilization  of  the  triad  tryptophan  cation  radical  in  APX  relative  to  CCP,  and  this  originates  only 
in  part  from  the  K+  cofactor  of  APX. 

2.  Methods 

The  x-ray  structure  coordinates  of  pea  cytosolic  APX  were  kindly  provided  by 
Professor  Thomas  L.  Poulos  [1 1, 12, 13];  PDB  accession  number  lAPX.  The  asymmetric  unit  of 
the  crystal  contains  four  identical  monomers,  which  are  labeled  A,  B,  C,  and  D  and  which  comprise 
two  homodimers.  These  were  utilized  separately  in  the  calculations.  The  crystal  structure  of 
recombinant  yeast  CCP  was  that  of  Goodin  and  McRee  [10];  PDB  accession  number  ICCA;  cf. 
Finzel,  Poulos,  and  Kraut  [15].  The  crystal  structure  of  M2301  CCP  (from  a  different  recombinant 
presentation)  was  that  of  Liu  et  al.  [16];  PDB  accession  number  2CEP. 

Electrostatic  calculations  of  indole  oxidation  were  performed  using  the  program  POLARIS  and 
the  Protein  Dipoles  Langevin  Dipoles  (PDLD)  method  [17-20].  The  PDLD  method  has  been 
applied  to  study  the  redox  properties  of  heme  in  cytochrome  c  [21-23],  of  photosynthetic  reaction 
centers  [24],  and  of  iron-sulfur  proteins  [25-28].  In  each  of  these  cases,  some  type  of  electrostatic 
representation  (partial  charges)  for  the  cofactors  within  the  proteins  were  derived  from  either  ab 
initio  or  semiempirical  methods. 
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The  PDLD  calculations  were  carried  out  using  an  IBM  RS6000  implementation  of  the  program 
POLARIS  [18].  The  redox  active  prosthetic  group  was  defined  to  be  region  I  and  was  in  this  case 
derived  from  a  3-methylindole  moiety.  The  protein  was  truncated  by  a  sphere  of  radius  rj  about  the 
centroid  of  region  I;  the  nonregion  I  protein  atoms  within  this  radius  define  region  n.  In  the  PDLD 
method,  solvent  water  was  modeled  microscopically  by  a  grid  of  orientable  but  positionally  fixed 
Langevin  dipoles.  The  Langevin  dipole  grid,  region  HI,  was  generated  within  a  sphere  of  radius  Tl 
centered  again  about  the  region  I  centroid.  In  all  cases,  rL  ^  rj.  The  Langevin  system  (region  III)  was 
further  radially  divided  into  an  inner  and  outer  region,  the  radius  of  the  inner  region  being  I2A.  The 
Langevin  dipole  grid  spacing  was  lA  and  3 A  in  the  iimer  and  outer  portions  of  region  DI, 
respectively.  Solvent  water  beyond  region  m  was  modeled  macroscopically  as  a  continuum 
dielectric  (region  IV)  using  a  Bom  expression  and  a  dielectric  constant  of  80. 

Starting  with  the  x-ray-derived  structures,  the  crystallographically  identified  atoms  assigned  as 
water,  O  atoms  were  deleted.  Hydrogen  atoms  are  added  using  standard  bond  lengths  and  angles 
[18].  Where  choices  of  H-atom  placement  were  not  obvious  (e.g. ,  a  side  chain  OH  group  or  histidine 
imidazole),  the  effects  of  varying  these  orientations  on  the  calculations  were  examined  and  minimum 
energy  positions  were  chosen.  For  protein  atoms,  heme,  and  K+,  standard  atomic  partial  charges, 
van  der  Waal’s  radii,  and  isotropic  polarizabilities  are  used  [18].  For  the  methylindole  cofactor, 
partial  charges  for  the  oxidized  and  reduced  states  for  region  I  atoms  were  derived  by  ab  initio 
methods  as  discussed  below.  Langevin  dipoles  were  deleted  at  grid  points  where  the  distance  to  any 
protein  atom  is  less  than  the  sum  of  the  atomic  van  der  Waal’s  radius  and  1.4A. 

The  PDLD  model  calculates  the  electrostatic  contribution,  AV,  of  the  protein  and  solvent 
surrounding  the  redox  active  prosthetic  group  to  the  free  energy  difference  of  the  oxidized  and 
reduced  states  of  the  redox  couple.  AV  is  therefore  a  solvation  energy  for  the  reduction.  AV  is  the 
sum  of  four  terms,  AVQp,  AVQa,  AVL,  and  AVB,  which  are,  respectively,  the  interaction  of  the 
region  I  atoms  with  the  fixed  partial  charges  of  the  protein  atoms,  the  induced  dipoles  of  the  protein 
atoms  arising  from  atomic  polarizabilities,  the  Langevin  dipole  grid,  and  the  continuum  dielectric. 
AVQa  was  produced  following  self-consistent  iteration  in  the  field  determined  by  protein  atoms. 
Solvent  water  was  modeled  by  orientable  dipoles  on  a  grid  defined  by  a  Langevin-Debye-type 
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equation.  Orientation  of  the  grid  dipoles  is  carried  out  iteratively  and  is  defined  by  the  electric  field 
resulting  from  protein  partial  charges,  cofactor  charges,  and  induced  protein  dipoles.  AVL  is  the 
energy  of  interaction  of  the  cofactor  with  the  field  defined  by  the  solvent  dipoles.  All  calculations 
unless  otherwise  specified  use  T2  =  22A  and  rL  =  25A. 

All  calculations  have  been  checked  for  convergence  with  respect  to  averaging  over  multiple 
Langevin  grid  positions,  and  iteration  to  achieve  the  self-consistent  protein  polarization  field  and 
Langevin  dipole  orientations.  For  CCP,  arginine-48  is  given  a  total  charge  of  +1  and  aspartate-235 
is  given  a  total  charge  of  -1.  For  APX,  the  analogous  arginine-38/aspartate-207  are  given  total 
charges  of  +1  and  -1,  respectively,  and  the  K+  atom  and  liganding  aspartate  187  are  given  total 
charges  of  + 1  and  - 1 ,  respectively.  All  other  residues,  the  heme,  and  charged  moieties  at  the  surface 
of  the  proteins  (heme  propionates,  residues  such  as  aspartate,  glutamate,  arginine,  and  lysine)  are 
modeled  as  neutral  [26]. 

The  free  energy,  AG,  for  the  oxidation  of  a  redox  cofactor  is  the  sum  of  the  intrinsic  energy  of 
oxidation  of  that  group  and  the  solvation  energy.  While  the  latter  is  approximated  by  the  PDLD 
method,  the  former  is  not  in  this  case  calculated,  and  therefore  only  differences  in  redox  potential 
due  to  differences  in  the  solvation  energy  (AAV)  are  modeled.  In  the  present  case,  given  the 
similarity  of  the  structural  environment  of  the  triad  tryptophan  in  CCP  and  APX,  AAV  is  likely  to 
result  from  the  relatively  subtle  differences  in  the  electrostatics  of  these  environments.  Redox 
potential  changes  can  be  calculated  for  tryptophan  in  APX  and  CCP  by  the  PDLD  method,  these 
being  obtained  by  using  AAV  for  the  two  proteins  and  the  conversion  factor  23.06  cal/mV. 
Molecular  dynamics  (MD)  were  also  applied  starting  with  the  x-ray  structures  and  using  the  program 
POLARIS  [18, 26, 27,  28]  in  order  to  provide  the  proper  averaging  over  configurations  generated 
in  both  redox  states.  Doing  so  provides  the  effect  of  protein  reorganization  and  yields  the  proper  free 
energy  within  the  linear  response  approximation  [18, 19].  MD  is  limited  to  atoms  lying  within  a 
sphere  of  radius  I2A,  centered  at  the  region  I  centroid.  In  addition,  during  MD,  Cartesian  constraints 
are  applied  in  APX  to  the  K+  and  aspartate- 187  06  atoms,  and  in  both  CCP  and  APX  at  the 
aspartate-235/207  06,  the  heme  liganding  histidine- 1 75/1 63  Ne,  and  tryptophan- 191/179  Ne  atoms. 
These  constraints  take  the  form  of  isotropic  atomic  harmonic  potentials  whose  minima  are  at  the 
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crystallographic  atomic  positions  and  whose  force  constants  (k)  are  50  kcal/A^  for  the  K+/ligating 
Oe  and  150  kcal/A^  for  the  tryptophan  Ne  aspartate  06  and  histidine  Ne  pairs.  MD  was  run  for  25  ps 
at  300  K.  PDLD  calculations  were  carried  out  on  50  MD  structures  taken  at  0.5-ps  intervals.  The 
results  were  averaged  together  with  a  t  =  0  value  to  5deld  MD-averaged  values  of  A  V. 

We  have  previously  reported  results  of  ab  initio  calculations  of  gas-phase  structures,  energies, 
harmonic  vibrational  frequencies,  infrared  and  Raman  intensities,  and  atomic  spin  densities  of 
3-methylindole  (skatole),  3-methylindole  neutral  radical,  and  3-methylindole  cation  radical  [9, 29]. 
These  calculations  employed  density  functional  theory,  the  Becke3LYP  functional  [30,  31]  as 
implemented  in  Gaussian  92/DFr  or  Gaussian  94  [32],  and  the  6-31G*  and  TZ2P  basis  sets  [33]. 
This  level  of  theory  has  been  shown  to  be  capable  of  providing  accurate  predictions  of  properties  for 
open  shell  molecules  [34, 35].  From  the  6-3 IG*  calculations,  the  atomic  partial  charges  used  for 
this  study  were  determined  using  the  electrostatic  potential  (ESP)  methods  developed  by  Kolhnan 
and  coworkers  [36,  37].  Specifically,  fitted  ESP  charges  were  obtained  with  the  Merz-Kollman 
option  in  Gaussian  94.  This  method  places  points  at  which  the  ESP  is  evaluated  on  four  shells 
around  each  atom  an  increasing  distance  from  each  nucleus  (1.2-,  1.4-,  1.8-,  and  2.0-times  the 
van  der  Waal’s  radius  of  that  nucleus).  The  points  are  placed  at  a  uniform  spacing  of  1  pt/Bohr  on 
each  surface.  The  atomic-centered  charges  are  fitted  to  the  values  of  ESP  with  a  Lagrange 
undetermined  multiplier  matrix  solution  that  permits  constraints  to  be  applied  to  the  system.  The 
only  constraint  used  here  was  that  the  sum  of  the  atomic  charges  equals  the  molecular  charge.  Since 
this  process  determines  the  atomic  partial  charges,  which  give  a  best  fit  to  the  ESP  on  a  surface 
surrounding  the  molecule,  they  are  expected  to  be  useful  when  imported  into  protein  electrostatic 
models,  such  as  is  implemented  in  the  PDLD  method.  In  order  to  assess  the  effects  that  the  polar 
protein  environment  might  have  on  the  ab  initio  properties  of  the  indole,  further  calculations  at  the 
Beck3LYP  6-3 IG*  geometry  were  carried  out  with  point  charges  placed  to  correspond  to  the  COj 
moiety  of  aspartate-235  of  CCP  (-0.75  for  each  oxygen,  +0.5  for  the  carbon)  and  for  the  amide 
groups  of  amino  acids  175  and  177  of  CCP  (+0.5  for  the  C=0  moiety,  ±0.45  for  the  NH  moiety). 
These  point  charge  positions  were  determined  by  fitting  the  3-methylindole  molecule  onto 
tryptophan-191  of  CCP  using  the  program  Xtalview  [38]  and  outputting  the  coordinates  of  the 
3-methylindole,  COj,  and  NHCO  moieties  for  input  into  Gaussian.  The  calculation  of  the  atomic 
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charges  in  the  presence  of  the  background  charge  distribution  was  done  using  the  Gaussian 
implementation  of  the  method  of  Hall  and  Smith  [39,  40]. 

3.  Results  and  Discussion 

Figure  1  shows  views  of  the  heme  active  site  and  tryptophan-aspartate-histidine  triad  containing 
regions  of  the  stmctures  of  CCP  and  APX.  Clearly,  the  structures  are  very  similar  to  one  another  in 
this  region,  despite  having  very  different  substrate  specificity.  The  similarity  of  structure  is  much 
less  surprising  in  view  of  the  ~33%  sequence  identity  between  the  two  peroxidases  [41].  Also 
shown  in  Figure  1  are  positions  of  Langevin  dipole  grid  points,  illustrating  the  similar  degree  of 
solvent  penetration  into  this  region  of  the  protein  for  both  enzymes.  Figure  2  shows  the  structure 
and  atom  numbering  for  the  analog  3-methylindole  cation  radical  used  to  model  the  tryptophan. 
Also  shown  are  the  ESP-derived  partial  charges  for  the  cation  radical  at  the  Becke3LYP/6-31G* 
level  of  theory.  Table  1  shows  partial  charges  for  3-methylindole  and  for  the  cation  radical  and  the 
values  translated  onto  the  triad  tryptophan  species  in  CCP  and  APX  (tryptophan-191  and 
tryptophan- 1 79,  respectively)  for  use  in  the  electrostatic  calculations.  Clearly,  the  cationic  character 
of  this  radical  is  distributed  over  the  whole  indole  moiety.  The  benzene  and  pyrrole  moieties,  for 
example,  have  total  charges  of  0.44  and  0.65,  respectively,  and  the  relative  charge  distribution  is  not 
altogether  different  from  unoxidized  3-methylindole.  Therefore,  the  analysis  of  electrostatic 
stabilization  by  either  the  nearby  triad  aspartate  in  CCP  and  APX  or  the  more  distant  K+  cation  in 
APX  is  more  complex  than  a  simple  Coulombic  interaction  between  two  point  charges.  In  APX 
(structure  A),  the  distances  between  the  K+  cation  and  the  trj^tophan  side  chain  atoms  range  from 
8.7A  (Cy)  to  1 1.5 A  (CC2).  Complicating  this  analysis  is  the  fact  that  the  K+  cation  is  more  properly 
viewed  as  an  aspartate- 187-K+  dipole,  with  an  overall  neutral  charge  (aspartate- 187  is  the  only 
charged  ligand  to  the  K+  binding  site).  This  dipole  will  have  a  less  straightforward  and  steeper 
distance  dependence  in  terms  of  interaction  with  a  cation  radical  at  the  triad  tryptophan. 

Table  2  shows  results  of  POLARIS  calculations  for  CCP  and  for  structures  A,  B,  C,  and  D  of 
APX.  Overall,  AV  for  CCP  is  larger  than  AV  for  any  of  the  APX  structures.  (Under  the  sign 
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Figure  1.  The  Active  Sites  of  Cytochrome-c-Peroxidase  (CCP)  and  Ascorbate  Peroxidase  (APX);  Structure  A.  Depicted  for 
Both  Proteins  Are  the  Distal  Histine,  Tryptophan  and  Arginine,  the  Heme,  and  the  Triad  of  Tryptophan,  Aspartate,  and 
Histidine.  Shown  for  APX  Is  Also  the  Aspartate-Potassium  Binding  Site.  For  Both  Proteins,  the  Locations  of  Langevin 
Dipoles  From  One  of  the  Dipole  Grids  Utilized  in  PDLD  Calculations  Are  Represented  as  Small  Spheres  (See  Section  2, 
Methods). 


(+0.18) 


(+0.17) 
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Figure  2.  Atom  Numbering  and  Becke3LYP/6-31G*  ESP  Charges  (in  Parentheses)  for 
3-Methylindole  Cation  Radical. 


radical  state  relative  to  the  unoxidized  state.)  The  range  of  AAV  values  (AV  for  APX  -  A  V  for  CCP) 
is  -3.5  to  -7.8  kcal/mole  for  the  four  available  APX  substructures.  This  translates  to  an  average 
redox  potential  difference  of  -279  mV  (range  -150  to  -339  mV).  The  primary  origin  of  this 
difference  is  in  the  first  Polaris  term,  AVQp,  which  is  between  16-20  kcal/mole  larger  in  CCP  than 
APX.  The  protein  polarization  term,  AVQa,  is  small  and  positive  in  APX,  and  for  CCP,  it  is 
-  6.6  kcal/mole.  The  value  AVL  also  contributes  significantly  to  AV  for  each  structure,  and  overall, 
the  sum  of  AVQa  and  AVL  serves  to  dampen  some  of  the  initial  difference  seen  in  AVQp.  AVB 
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Table  2.  PDLD  Results  With  X-ray-Derived  Structures’ 
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is  relatively  invariable  for  each  structure.  Details  of  the  origin  of  these  terms  in  CCP  and  APX 
follow. 

An  important  observation  of  this  study  is  the  variation  in  the  electrostatic  energies  for  the  four 
examples  of  the  APX  structure  within  the  crystal  asymmetric  unit  (the  overall  root  mean  square 
positional  deviation  for  Ca  atoms  is  ^  0.4A  for  any  given  pairwise  comparison  of  the  four  structures) 
[11-13].  Differences  in  these  four  structures  most  likely  result  from  packing  effects  and  may  reflect 
average  snapshots  of  the  range  of  conformers  available  to  the  protein  in  solution.  In  other  studies 
using  the  PDLD  method  to  study  redox  potentials,  great  benefit  in  terms  of  method  accuracy  was 
realized  from  generating  a  family  of  structures  by  MD  from  any  x-ray  derived  structure  [18, 1 9, 20, 
26, 27, 28].  The  improvement  might  in  any  given  case  result  from  adjusting  structural  elements  that 
are  not  energetically  reasonable,  sampling  a  range  of  thermally  accessible  structures  rather  than  the 
x-ray-derived  average,  improving  the  convergence  of  the  iterative  aspects  of  the  PDLD  model 
(especially  in  the  solvent  mode),  sampling  multiple  orientations  for  the  placement  of  freely 
orientable  hydrogens,  inclusion  of  protein  reorganization  energy,  or  all  of  the  above.  Figure  3  shows 
AV  determinations  for  51  structures  generated  during  MD  for  CCP  and  APX  (structure  A,  t=0  and 
25  ps  of  dynamics  with  0.5-ps  snapshots;  structures  B,  C,  and  D  had  similar  profiles).  Figure  4 
shows  the  structures  of  the  snapshot  samples  for  CCP.  MD  averaged  values  for  all  of  the  structures 
are  shown  in  Table  3.  In  each  case,  AV  for  a  given  structure  decreased  upon  MD  averaging.  In 
addition,  the  spread  in  AV  values  for  the  four  APX  structures  decreases,  and  AAV  for  the  average 
of  all  four  APX  structures  vs.  CCP  decreases  to  -7.6  kcal/mole  or  -330  mV. 

The  fully  averaged  (MD,  multiple  structures)  values  of  AAV  provide  for  the  prediction  that 
tryptophan-179  in  APX  is  330  mV  more  difficult  to  oxidize  to  the  radical  cation  than  tryptophan-191 
in  CCP.  Mondal,  Fuller,  and  Armstrong  [42]  have  reported  that  the  midpoint  potential  of  this 
oxidation  in  CCP  is  +740  mV  vs.  S.H.E.  This  places  the  estimated  potential  for  APX  at  greater  than 
1  V,  which  is  higher  than  the  midpoint  potential  for  indole  water  [43].  We  cannot  at  this  time  rule 
out  all  other  options  for  the  absence  of  an  observed  tryptophan  radical  cation  in  APX.  The  radical 
could  exist  but  have  a  lifetime  too  short  to  observe  under  the  experimental  conditions  employed  to 
date,  the  porphyrin  n-cation  radical  that  is  observed  could  simply  be  unusually  stable  in  APX  relative 
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Figure  3.  Values  of  AV  for  APX  (Structure  A)  and  CCP  Over  25  ps  of  Molecular  Dynamics. 

to  CCP,  or  there  could  be  a  kinetic  barrier  to  oxidation  of  tryptophan  in  APX  sufficient  to  divert  the 
second  oxidizing  equivalent  elsewhere.  Also,  the  minimum  value  of  the  midpoint  potential  shift 
from  CCP  to  APX  sufficient  to  shunt  the  second  oxidizing  equivalent  of  HjOj  to  another  site  is  not 
known.  However,  these  calculations  do  show  that  electrostatic  differences  alone  between  CCP  and 
APX  are  large  enough  to  explain  the  thermodynamic  instability  of  this  radical  species  in  APX  on 
reaction  with  H2O2.  This  conclusion  is  consistent  with  the  observation  that  cation  binding  sites  can 
be  created  in  the  active  site  of  CCP  by  cavity-forming  site-directed  mutants  [44-53],  Given  a 
tryptophan  site  electrostatically  predisposed  for  stabilizing  a  cation,  it  is  also  not  surprising  to  note 
that  the  structure  of  CCP  is  little  changed  upon  oxidation  to  compound  ES  [54]. 
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Figure  4.  Fifty-One  Superposed  ‘‘Snapshot”  Structures  From  a  25-ps  Molecular  Dynamics  Run  for  CCP.  PDLD  Calculations  for 
These  Structures  Are  Utilized  to  Determine  Structure-Averaged  Results. 


Table  3.  PDLD  Results  From  MD  Averaging  of  X-ray  Structures 
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To  trace  the  structural  origin  of  the  electrostatic  differences  seen  in  the  calculations  of  AAV,  we 
examined  the  individual  protein  residue  contributions  to  AVQp,  since  this  term  dominates  the 
difference  in  AAV.  Figure  5  displays  a  profile  of  per  residue  contributions  to  AVQp  for  both  CCP 
and  APX  (structure  A  is  shown;  structures  B,  C,  and  D  are  all  very  similar).  The  homology  between 
the  two  proteins  is  obvious.  Prominent  features  (moving  linearly  in  sequence  space),  which 
contribute  about  equally  in  CCP  and  APX,  are  from  the  distal  arginine  cation  (48  in  CCP,  38  in 
APX),  the  trio  of  backbone  carbonyl  groups  firom  residues  172, 175,  and  177  in  CCP  and  160, 163, 
and  165  in  APX,  and  the  active  site  aspartate  (235  in  CCP,  207  in  APX).  The  three  major 
differences  between  CCP  and  APX  are  the  larger  contribution  of  the  1 72/ 1 75/1 77  backbone  carbonyl 
trio  in  CCP  relative  to  the  1 60/ 1 63/ 1 65  trio  in  APX  (18.7vs.  13.2  kcal/mole),  the  significantly  larger 
contribution  attributed  to  methionine  230  in  CCP  relative  to  leucine  203  in  APX,  and  the  presence 
in  APX  of  the  aspartate  187/K+  dipole  with  no  equivalent  in  CCP.  Table  4  lists  the  values  AVQpi 
for  each  of  these  interactions  as  well  as  other  smaller  contributions. 

All  the  significant  contributions  noted  here  are  within  15A  of  the  triad  tryptophan.  Calculations 
run  with  rj  set  at  (rL  =  25A)  give  AAV  values  essentially  identical  to  those  reported  with  r2  =  22A 
(data  not  shown),  confirming  the  relatively  local  nature  of  the  origin  of  AAV.  The  distal  arginine 
(48  in  CCP,  38  in  APX)  resides  in  a  partially  solvent-filled  cavity  on  the  distal  side  of  the  heme  with 
a  nearest  approach  of  ~lk  to  the  relevant  tryptophan.  Calculations  with  this  residue  modeled  as 
electroneutral  yield  AAV  values  essentially  identical  to  those  which  use  a  charge  of  +1  on  this 
residue  (data  not  shown).  Methionine  230  contributes  to  AVQp  from  both  the  polar  side  chain  (Sj 
is  ~Sk  from  tryptophan  191)  and  from  the  backbone  carbonyl  (which  oxygen  is  ~2>k  from  the 
tryptophan).  Leucine  203  in  APX  has  a  nonpolar  side  chain  and,  in  addition,  has  a  less  favorable 
orientation  of  the  corresponding  backbone  carbonyl  relative  to  CCP.  We  have  performed 
calculations  using  the  structure  for  the  CCP  mutant  where  methionine  230  is  replaced  by  isoleucine 
(M230I)  [16].  Results  for  the  x-ray  structure  and  for  MD  average  structures  are  presented  in 
Tables  2  and  3.  Using  the  x-ray  structure,  there  is  a  8  kcal/mole  decrease  in  AVQp  in  M230I-CCP 
relative  to  CCP,  which  contributes  to  a  AAV  of  -  3. 1  kcal/mole  or  - 1 36  mV.  MD  averaging  for  this 
stmcture,  however,  reduces  AAV,  now  positive,  to  only  0.3  kcal/mole  or  13  mV.  Experimentally, 
although  the  kinetics  for  reaction  of  H2O2  were  altered  by  the  M230I  mutation  [16],  a  tryptophan 
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Figure  5.  AVQiii  (the  Contribution  to  AVQp  for  the  Redox  Transition  From  Each  Amino  Acid  i).  The  Last  Entries  for  CCP  and 
APX  (Monomer  A)  Are  the  Heme  and  the  Heme  and  K-)-  Atom,  Respectively.  Note  That  for  the  Sign  Convention 
Employed  Here,  a  Positive  Value  of  a  AVQpi  Is  Electrostatically  Stabilizing  to  a  Cation  Radical  at  the  Triad 
Tryptophan. 


Table  4.  Contributions  to  AYQu 
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radical  is  reported  to  form  in  this  mutant.  So  while  methionine  230  provides  stabilization  to  the 
tryptophan  radical  in  CCP  relative  to  leucine  203  in  APX,  mutation  at  this  site  results  in  relatively 
little  change  in  overall  electrostatic  stabilization  of  a  tryptophan  radical.  This  illustrates  the  point 
that  a  given  moiety  (e.g.,  an  amino  acid  side  chin)  can  contribute  to  the  electrostatic  destabilization 
of  a  redox  process,  but  upon  mutation,  the  system  (protein  structure,  polarity,  and  solvent)  can 
“relax”  to  (partially)  accommodate  the  change. 

The  most  striking  factor  noted  in  the  structure  of  APX  relative  to  CCP  was  the  presence  of  the 
potassium  site.  This  cation  site  has  been  proposed  [11, 12, 13,  55]  as  a  possible  mechanism  for 
electrostatic  destabilization  of  a  tryptophan- 179  radical  in  APX.  Consistent  with  this  proposal,  the 
electrostatic  contributions  in  Table  4  of  the  K+  and  the  ligand  aspartate- 1 87  are  significant  and  find 
no  analog  in  CCP.  However,  the  ion-pair  forms  a  dipole  (see  Table  4),  which,  while  favorably 
oriented  to  destabilize  a  cation  at  tryptophan- 1 79,  is  some  9- 1 OA  away,  with  both  intervening  solvent 
and  polypeptide  medium  in  between  (Figure  1)  [11-13].  The  PDLD  model  can  evaluate  this 
interaction  without  assumption  of  a  bulk  dielectric  constant  for  this  composite  medium  [56] ,  and  we 
have  examined  the  effects  of  eliminating  the  charges  on  the  aspartate-K+  pair  in  APX.  Results  for 
these  calculations  are  reported  in  Table  2  for  x-ray-derived  stractures  and  in  Table  3  for 
MD-averaged  calculations.  For  the  x-ray-derived  structures  without  MD  averaging,  the  data  in 
Table  2  show  that  the  aspartate-K+  ion-pair  contributes  relatively  little  to  destabilizing  the 
tryptophan  cation  radical.  Without  the  charged  ion-pair,  the  average  AAV  was  -260  mV,  only 
19  mV  smaller  than  the  corresponding  calculation  with  the  aspartate-K+  dipole  present.  As  before, 
a  spread  of  values  is  observed  for  structures  A-D,  and  MD  averaging  again  reduces  the  spread 
between  these  structures.  After  including  MD  averaging  (MD,  stractures  A-D),  we  obtained  AAV 
between  the  ion-pair  charge  deleted  structures  and  CCP  of  -7.1  kcal/mole,  which  corresponds  to 
-308  mV,  only  ~20  mV  different  from  AAV,  which  includes  the  ion  pair.  From  ttiis  result,  it  can 
be  concluded  that  the  aspartate-K+  pair  does  play  a  role  in  destabilizing  a  cation  radical  at 
tryptophan- 179  in  APX  but  is  only  one  of  several  factors.  Of  significant  interest,  a  recent  study 
reported  the  engineering  of  a  K+  binding  site  m  CCP  by  site-directed  mutagenesis  at  the  analogous 
position  to  that  of  APX  [13]  resulting  in  a  significant  destabilization  of  the  tryptophan- 191  radical 
[57-60]. 


19 


Absolute  errors  in  the  calculations  presented  here  are  expected  to  be  significant,  but  are  reduced 
substantially  by  the  problem  posed  (comparing  electrostatic  stabilization  of  a  single  species  in  two 
very  similar  environments).  Modeling  of  the  aspartate-tryptophan  interaction  is  at  this  point  purely 
electrostatic,  without  consideration  of  hydrogen  bonding.  The  heme  is  treated  as  overall 
eletroneutral,  but  of  course  in  CCP  this  group  also  undergoes  oxidation  from  a  formally  cationic  to 
a  neutral  state.  Molecular  dynamics  required  unsophisticated  Cartesian  constraints  on  the  histidine 
heme  ligand,  triad  aspartate-tryptophan,  as  well  as  on  the  aspartate-K+  pair.  With  the  exception  of 
the  latter,  all  of  these  approximations  are  common  to  both  stmctures  and  therefore  are  not  expected 
to  translate  into  significant  errors  in  AAV.  While  the  ab  initio  methods  used  to  obtain  ESP  charges 
are  quite  sophisticated,  they  represent  an  in  vacuo  state  of  the  3-methylindole  molecule.  Examining 
the  effect  of  the  highly  charged  CCP  environment  (specifically,  the  charge  on  aspartate  235  and  the 
amide  groups  from  residues  175  and  177)  on  these  charges  revealed  only  very  small  changes  in 
charge  distribution  (Table  1).  Calculations  of  AV  and  AAV  were  similarly  affected  very  little  by 
these  different  charges  (data not  shown).  Although  the  solvent  model  employed  in  the  determination 
of  AVL  is  a  microscopic  model,  it  falls  short  of  a  full  modeling  of  the  aqueous  environment. 
Figure  6  shows  the  number  distribution  and  energy  distributions  of  individual  dipoles  contributing 
to  AVL.  Clearly,  modeling  of  solvent  in  a  reasonably  sophisticated  manner  will  improve  the 
accuracy  of  any  calculation  of  AAV. 

Calculations  of  the  electrostatic  environment  of  tryptophan- 179  of  APX  relative  to 
tryptophan- 1 9 1  of  CCP  have  been  reported  previously  [11-13].  A  continuum  model  was  used  with 
dielectric  constants  of  either  2  or  9.  Charge  parameters  for  the  tryptophan  cation  radical  were  not 
specified.  These  studies  clearly  demonstrated  a  more  positive  ESP  in  APX  relative  to  CCP  and  also 
suggested  a  partial  role  for  the  K+  atom  in  establishing  this  difference  in  potential.  However, 
quantitative  interpretation  of  these  results  must  suffer  from  the  arbitrary  choice  of  dielectric,  which 
can  vary  enormously  for  different  interactions,  with  or  without  significant  presence  of  solvent  water 
(see  discussion  in  Jensen,  Warshel,  and  Stephens  [26]  and  Warshel,  Papzyan,  and  Muegge  [56]). 
Even  relative  comparisons  (e.g.,  AAV  for  APX  vs.  CCP)  may  be  difficult  to  derive  using  continuum 
methods,  since  the  specific  between  group  effective  dielectric  constants  [61]  may  be  different  for 
different  interactions  (e.g.,  triad  aspartate-tryptophan  vs.  K+-tryptophan).  The  present  calculations 
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Figure  6.  Upper  Plot:  Cumulative  Occurrence  of  Langevin  Dipoles  Moving  Radially  Away 
From  the  Geometric  Centroid  of  the  Triad  Tryptophan  in  CCP  and  APX 
(Structure  A).  Lower  Plot:  Energetic  Contribution  to  AVL  for  CCP  and  APX  of 
Each  SA  Shell  (O-SA,  3-6A,  etc.). 
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expand  on  the  previous  work  in  part  by  including  realistic  modeling  of  the  tryptophan  radical,  protein 
polarization,  MD  structure  averaging,  and  a  microscopic  solvent  model,  and  provide  quantitative 
estimates  of  the  difference  in  redox  potential  for  the  active  site  tryptophan  between  CCP  and  APX, 
as  well  as  the  magnitude  of  the  principle  factors  that  contribute  to  this  difference. 
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